Realization of quantum process tomography in NMR 
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Quantum process tomography is a procedure by which the unknown dynamical evolution of an 
open quantum system can be fully experimentally characterized. We demonstrate explicitly how 
this procedure can be implemented with a nuclear magnetic resonance quantum computer. This 
allows us to measure the fidelity of a controlled-NOT logic gate and to experimentally investigate 
the error model for our computer. Based on the latter analysis, we test an important assumption 
underlying nearly all models of quantum error correction, the independence of errors on different 
qubits. 



PACS: 03.67.-a, 06.20.Dk, 76.60.-k 

I. INTRODUCTION 

Experimental characterization of the dynamical behav- 
ior of open quantum systems has traditionally revolved 
around semiclassical concepts such as couplingstrengths, 
relaxation rates, and phase coherence times How- 
ever, quantum information theory tells us that so few pa- 
rameters barely even begin to represent the full dynamics 
which quantum systems are capable of; for example, the 
evolution between two fixed times of a two-level quantum 
system (a qubit), coupled to an arbitrary reservoir, is de- 
scribed by twelve real parameters ||. For a two qubit 
system, this number grows to 240, and in general, for n 
qubits it is 16" - 4". 

Of course, not all of these parameters are relevant for 
every physical process, and understanding the physics 
can boil down this huge parameter space to just a few 
important numbers. On the other hand, when the physi- 
cal origins of a quantum process are not understood, or in 
question, it is invaluable to know that all these param- 
eters can, in principle, be experimentally measured - 
albeit with exponential effort (in n) — using a procedure 
known as quantum process tomography. This method is 
a direct quantum extension of the classical concept of 
"system identification" which enables control of dynam- 
ical systems, and provides powerful measurement tech- 
niques for characterizing unknown systems. For closed 
systems, process tomography is the analog of determin- 
ing the truth table of an unknown function; quantum- 
mechanically, this means the unitary transform involved 
in a quantum computation. 

The procedure for quantum process tomography has 
been described in detail in the literature ^-^]. It re- 
lies upon the ability to prepare a complete set of quan- 
tum states pj as input to the unknown process St, and 
the ability to measure the density matrices of the out- 
put quantum states from the process £t{p)- As a result, 
one obtains a complete "black box" characterization of 



St- By systematically repeating this procedure for in- 
creasing process time t, one can also obtain a complete 
master equation description of the process |6) . 

This procedure seems to provide an ideal way to char- 
acterize and test quantum gates on a small number of 
qubits — for example, those realized using nuclear mag- 
netic resonance (NMR) techniques However, to 
date, this has not been practiced because of two prob- 
lems: first, it is desirable to explicitly include the max- 
imally mixed state as one of the Pj's, so as to directly 
subtract its contribution, but such a state cannot be pre- 
pared by unitary actions alone. Second, one can only 
measure traceless observables, and thus one cannot ob- 
tain the entire output density matrix of a process. 

Here, we show explicitly how these hurdles can be over- 
come, and demonstrate the use of these methods in ad- 
dressing two outstanding issues in quantum computation: 
the fidelity of a quantum logic gate, the controlled-NOT 
gate, and the validity of the independent error model 
which is widely assumed for quantum error correction 
and fault tolerant computation. 

The paper is divided into two major sections, which 
describe the theory, and our experiment. We begin with 
a review of quantum process tomography (QPT), then we 
explain our extensions to the basic procedure which en- 
able QPT with NMR, and present theoretical models for 
decoherence in NMR. We then describe our experiment 
and procedure, and our two main experimental results. 

II. THEORY 

A. Quantum process tomography 

Quantum process tomography, as introduced by 
Chuang and Nielsen j|, can be summarized as follows. 
A general quantum operation on a quantum state p is 
a superoperator £, a linear, trace-preserving, completely 
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positive map We consider only the case where 

£{p) has the same dimension as p. A common form for 
£ is the operator-sum representation, 



£{p) = Y J A k pA\, 



(1) 



where J^k ^i^-fc = ® ne can easily demonstrate the 
transformation to a fixed-basis expansion of the quan- 
tum operation such that the information about the pro- 
cess lies in a set of coefficients Xmn 

instead of a set of 

operators. In this form, £ maps an initial density matrix 
p according to 



£(P) = ^XmnAmpAl 



In this expression, the A m are a basis for operators on 
the space of density matrices, and x is a positive Her- 
mitian matrix with elements Xmn- In the following, we 
will omit explicit summations and adopt the convention 
that repeated indices are summed over. Because of the 
restriction that £ must preserve the trace of p, \ contains 
TV 4 — iV 2 independent parameters for an N dimensional 
system (for n spins, N = 2"). 

Let the N 2 matrices pj be a basis for density matrices. 
Applying £ to this basis gives 



£{Pj) — ^jkPk 



(3) 



Using quantum state tomography |p.2| , one can experi- 
mentally determine Xjk, which fully specifies £. To de- 
termine x from A, we define (3 by 



AmPjA n — (3jk n pk , 



(4) 



which is fully specified by the choice of A m and pj . Then 
one can show that 



Mjk Xmn 



A 



jk ■ 



(5) 



We may think of (3 as a matrix and A and x as vectors, 
with mn a composite column index and jk a compos- 
ite row index; then fix = A. Using the pseudoinverse k 
of (3 (also called the Moore-Penrose generalized inverse) , 
which may be computed independent of measurement re- 
sults, we have 

X = k\. (6) 



B. QPT in NMR 

Implementing QPT in high temperature NMR presents 
two main obstacles. First, the part of the system which 
can be manipulated by unitary operations is but a small 
deviation from a maximally mixed state. Second, one 
can only measure traceless observables, and the scaling 



of these observables with respect to the traceful part of 
the system is not immediately known. Fortunately, both 
of these hurdles can be overcome. 

The Hamiltonian for solution NMR is well approxi- 
mated by 



H 



hj. 



jk- 



ZjZk 



where 



(7) 



(8) 



is the Pauli operator for the jth spin along the direction 
(2) of the magnetic field, u>j is its Larmor frequency, and the 



sign of the first term is chosen so that the ground state 
is spin up. The second term represents a small scalar 
coupling between spins known as the J-coupling. 

The initial state of an NMR quantum computer is the 
thermal equilibrium state with density matrix 



e -H/k B T 
tr ( e -if/fc B T) 



= cl + A c 



(9) 



where c = 2~ n for a system of n spins (so that tr(poo) = 
1) and Aqo is traceless. In general, we refer to the trace- 
less part of any density matrix p as the deviation density 
matrix. If we neglect the small coupling between spins 
and consider high temperature with respect to the Lar- 
mor frequencies, the equilibrium density matrix is ap- 
proximately 



2k B T 



(10) 



By performing unitary operations on the system, we 
may rotate the ujjZj deviation, allowing us to prepare 
a linearly independent set of inputs pj. Because all of 
the interesting information will lie in this deviation, we 
would like to be able to prepare an input with no devi- 
ation (i.e., the maximally mixed state), so that we may 
directly subtract the contribution of the identity. Such an 
input clearly cannot be prepared by unitary means. How- 
ever, the maximally mixed state can be prepared from the 
thermal state by applying a ^ pulse followed by a uni- 
formly spatially varying RF pulse. Different parts of the 
sample experience different amounts of rotation, so that 
the ensemble is effectively dephased to the maximally 
mixed state. 

The second difficulty can be handled by calibrating 
measurements to the known initial state of the system. 
By state tomography, we can measure A = dAoo, where 
the constant d is arbitrary — it depends in detail on 
the gain of the amplifiers, the efficiency of the RF coils, 
etc. But the initial deviation density matrix is known 
to be c- 



2k B T 



Hence we may calibrate our apparatus by 

determining the value d that makes the measured A/d 
closest to the theoretical Aqq. In practice, we effect this 
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calibration by normalizing all measurements (integrals 
of spectral peaks) for spin j to a reference measurement 
taken on that spin (a ^ pulse applied to the thermal 

state) multiplied by c ■ 



C. Decoherence primitives 

Decoherence of a single NMR spin in solution is gen- 
erally characterized by two well-known processes, am- 
plitude damping and phase damping. Here, we present 
these processes, as well as generalizations to multiple 
spins. To facilitate the description of multiple simultane- 
ous processes, we describe amplitude damping and phase 
damping as semigroups in terms of their generators. 

Assuming that a quantum operation can be de- 
scribed by a norm continuous one-parameter semigroup, 
given a superoperator £ t satisfyingthe stationarity and 
Markovity condition £ s £t = £ s +t ]lq ], its generator Z is 
defined as 111 



Zip) = lim 

t^o 



£«0) - p 



t 



(ii) 



The generators have the advantage that they easily de- 
scribe the simultaneous action of multiple processes. In 
particular, we have the Trotter formula |l4|Jlq|, 



£, 



(1+2) 



= lim 



\°t/n C t/n) ' 



where 



Z (l+2) = z (l) + Z {2) 



(12) 



(13) 



with obvious generalization to more than two processes. 

Furthermore, we may return to the superoperator form 
by exponentiating the generator pjj: 



£ = lim ( 1 - -Z 

zt 



(14) 
(15) 



This exponentiation may be defined by its Taylor series 
in the case of interest where Z is bounded. To imple- 
ment exponentiation in practice, it is useful to work with 
a manifestly linear representation, i.e., that of Eq. (^). In 
this case, a superoperator and its generator can be rep- 
resented as matrices which have the property that com- 
position of superoperators corresponds to matrix multi- 
plication. Then the procedure of Eq. (|l4|) corresponds to 
matrix exponentiation. 

We now consider the single-qubit versions of phase and 
amplitude damping. For further discussion of these pro- 
cesses, including operator-sum representations, sec [ |16| . 
We write the density matrix of a qubit as 



Poo Poi 
Pio Pn 



(16) 



Phase damping can be thought of as a consequence of 
random phase kicks. It acts on a qubit as 



trip) 



Poo 

" 7 'pio 



e 7 *poi 
Pn 



(17) 



for some damping rate 7; i.e., it has a generator which 
acts as 



{ ° 

\Pio 



Poi 




(18) 



Ordering the matrix elements of p in the vector 
(poo; Poii P10, Pn) T (which effectively defines the basis 
used in Eq. (||)), we have 



rPD 
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(19) 



Generalized amplitude damping is a process whereby 
the qubit may exchange energy with a reservoir at some 
fixed temperature. It acts on a qubit as 



-GAD 



(P) 



fclPOO + k 2 Pll 

„-rt/2 „ 
e ' pio 

-rt\ 



e -rt/2 

hPOO 



P01 

&4P11 



(20) 



where = (1 — n)(l - e _lt ), fc 3 = n(l - e~ 14 ), k\ = 
1 — &3, and &4 = 1 — k 2 , with T a damping rate and n a 
temperature parameter. Its generator acts as 



Z 



GAD 



= -r 



I n n-Y 
i 
i 



V 



(21) 



-n 



1 



Computing the fixed point of this process, we find that 
we may interpret it as occuring at a temperature 



k B T = 



AE 
logi^ 



(22) 



where AE is the splitting between the ground and excited 
states of the system. 

Heuristically, we expect phase damping to describe the 
T2 (dephasing) process of NMR. Similarly, generalized 
amplitude damping should describe the T\ process of 
NMR, whereby a spin relaxes to align with the magnetic 
field. Note from the geometry of the Bloch sphere that 
amplitude damping necessarily includes loss of phase in- 
formation, so it also contributes to the Ti process. How- 
ever, in typical systems, T\ is much longer than T2, so 
the contribution of amplitude damping to dephasing is 
small. 

The simplest extension of these processes to a multiple- 
spin system is to allow the individual processes to act 
independently on each qubit. However, this is far from 
the most general scenario: the decoherence might be cor- 
related. Although we have been unable to find a suitable 
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correlated generalization of amplitude damping, we have 
incorporated a simple model of correlated phase damping 
due to Zhou into our analysis ]l7j] . In this model, each 
spin experiences a random phase shift with zero mean. 
The amount of damping and the extent of correlation is 
defined by the covariance matrix of the phase shifts. By 
explicit calculation, one may easily find that this leads 
to a generator of the form 



rCPD 



= diag[0, -72, -71, -(71 + 72+73), 
-72,0, -(71 +72 -73), ~7i, 
-7i, -(7i +72 -73), 0,-72, 
-(71 +72 +73), -71, -72,0], 



(23) 



where 71 and 72 may be interpreted as rates for inde- 
pendent phase damping on spins 1 and 2, and 73 may 
be interpreted as a rate for correlated damping. Indeed, 
in the case 73 = 0, this model reproduces independent 
phase damping. 

Our completed model has the generator 



Z = z- 



jjCPD _|_ ^jGAD 



z. 



GAD 



(24) 



where a subscript indicates which spin is acted on and 
Z J is the generator corresponding to the Hamiltonian 



hJ 



(25) 



In this Hamiltonian, the Zeeman terms are dropped from 
Eq. (Q) because we work in a frame rotating at the Lar- 
mor frequencies of the two spins. To produce the model 
superoperator at any given time, we simply exponentiate 
Eq. (g). 



III. EXPERIMENT 



be captured by many distance measures 
the minimum gate fidelity, defined as 



F 



mm 

IV>> 



— for example, 
(iP\tf£(\i>)(i>\)U\iP). (26) 



F represents the smallest possible overlap between a state 
acted on by U and the same state acted on by £. Note 
that minimizing over pure states is sufficient because the 
fidelity is convex and because an arbitrary density matrix 
can be written as a convex sum of pure states. 

Using QPT, we measured the minimum gate fidelity 
for a controlled-NOT operation, a well-known member of 
universal gate sets, with 



(\ 0" 

10 

1 

\o 1 0, 



(27) 



The NMR pulse sequence chosen to implement this gate 
was 



n a-, >-,}•:../ 1 -L ) jo .y,. 



(28) 



where time goes right to left, Xj denotes pulse 
on spin j, a bar denotes an inverse pulse, d(-) denotes a 
time evolution during which no pulses are applied, and 
J = 215 Hz is the coupling constant for the scalar J- 
coupling between the spins. Procedures similar to those 
used in fl8|| were adopted to implement this pulse se- 
quence. The pulse lengths were approximately 10 — 20 /is, 
sufficiently fast that the effect of J-coupling could be ne- 
glected during a pulse. Using quantum process tomogra- 
phy, we determined x f° r t ms process, effectively deter- 
mining £. This result is shown in Fig. [l]. 



We have implemented QPT on a system of two spins 
in an NMR apparatus. The experiments were performed 
at the IBM Almaden Research Center using an Oxford 
Instruments wide-bore magnet and a 500 MHz Varian 
Unity Inova spectrometer with a Nalorac triple resonance 
HFX probe. The sample was approximately 0.5 mL of 
200 millimolar isotopically labeled chloroform ( 13 CHCl3) 
in d6-acetone, where the two qubits were the proton (first 
spin) and carbon (second spin) . This volume was chosen 
to make the sample fairly short, as this gave the most 
effective gradient pulse characteristics. 



A. Results: Gate fidelities 

One application of quantum process tomography is to 
diagnose the accuracy of a quantum computation. When 
one applies a series of gates to perform the unitary oper- 
ation U, the computer will actually perform a superoper- 
ator £, which one hopes is close to U. That closeness can 





FIG. 1. x matrices for the controlled-NOT gate. The ma- 
trix on the left was experimentally measured by QPT, and 
that on the right is theoretical. Only the magnitudes of the 
matrix elements are shown; the phases also corresponded well 
between theory and experiment. The matrix element in the 
far left corner connects II <-> II, and the elements are ordered 
as II,XI,...,YZ, ZZ. 

Given Xi we ma Y calculate gate fidelities. For the 
controlled-NOT gate, numerically minimizing over input 
states gives -Fcnot = 

0.80 + 0.04. 
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This figure may be thought of as a rough benchmark 
for the quality of gates that can be implemented us- 
ing current NMR quantum computers. However, sev- 
eral caveats should be addressed. Foremost, the prepa- 
ration and readout steps performed to do the process 
tomography contribute to much of the error. Perform- 
ing process tomography on a null computation (U = I), 
we find Fi — 0.90 ± 0.03. The primary sources of error 
are most likely imperfect state preparation and imperfect 
state tomography due to imperfect pulse calibration and 
inhomogeneity of the RF fields JTs) ] . 

Also, in a long computation, experimental results sug- 
gest that there may be a significant cancellation of er- 
rors |Q. Typically, the major contribution to the error 
introduced by an individual gate will be largely due to 
systematic errors rather than fundamentally irreversible 
decoherence. If the pulse sequence exhibits some degree 
of symmetry, these systematic errors may at least par- 
tially cancel. Thus the fidelity of a long sequence of 
gates cannot be deduced from the individual gate fideli- 
ties, though one would certainly expect it to be subaddi- 
tive. 

Of course, the minimum gate fidelity is a pessimistic 
standard; it may be more reasonable to consider the av- 
erage gate fidelity. To numerically calculate the aver- 
age fidelity, we must be able to sample uniformly from 
quantum states according to an appropriate probability 
measure. Although this problem is subtle for the case of 
mixed states, there is a straightforward choice for pure 
states: we take the unitary transform of a fixed state, 
where the unitary operator is chosen according to the 
Haar measure (the unique invariant measure on a Lie 
group). Using the procedure described in ^0|, we cal- 
culated average fidelities with 10 5 randomly chosen uni- 
taries, giving F c a ^ T = 0.955 and Ff g = 0.960. Thus, 
the average fidelity is significantly better than the worst 
case; the error rate 1 — F improves by about one order 
of magnitude. 

B. Results: Decoherence characterization 

To characterize the decoherence occuring in the chlo- 
roform system, we used QPT to measure \ over the var- 
ious relevant time scales. For an NMR spin, the rele- 
vant scales are jj , T\ , and T2 |2^| . We thus sampled at 
times given by (^Y jj for integers j E [0,4], (^pTi for 
j E [-1,4], and (|) J T 2 for j £ [0,9], where J = 215 Hz 
is the known value of the J-coupling, T\ — 20 s is a time 
scale on the order of T\, and T2 = 0.5 s is an intermediate 
time scale on the order of T2 . 

By executing process tomography, we are able to de- 
termine x as a function of time. However, x is a large, 
complex collection of numbers that cannot be easily in- 
terpreted. To better understand the results, we fit the 
data to a model process, Xm, which hopefully provides 



a reasonable description of the relevant physics. Specifi- 
cally, we use the model which results from exponentiating 
Eq. (§]). 

For each x, we determined the closest fit to our model 
by numerically minimizing tr[(x m — xV(Xm — x]L where 
Xm is a model superoperator derived from Eq. ( |24| ) . This 
figure of merit was on the order of 5% for most fits. The 
experimentally measured x matrices as well as the cor- 
responding fits Xm are shown in Fig. for three delay 
times. 

Fitting the various rate parameters of Xm as a function 
of time allows us to characterize the process. These fits 
are shown in Figs. ^-j^. The error bars in these plots are 
computed numerically, and are derived solely from the 
statistics of the measurements and the fitting procedure. 




FIG. 2. x matrices for decoherence at various times. Ma- 
trices on the left are experimentally measured using QPT, 
and those on the right are fits to Eq. (p^). Only the mag- 
nitudes of the matrix elements are shown, but their phases 
were also described well by the model. From top to bot- 
tom, t = 0.065 s, 0.5 s, 20 s. At t = 0.065 s, the time was 
chosen to be an integer number of J-coupling periods, which 
explains the lack of IZ, ZI, and ZZ terms in the corners. At 
t = 0.5 s, note that the IZ <-> ZI terms arise from a combina- 
tion of J-coupling and independent phase damping, not from 
correlated phase damping. At t — 20 s, note that the nonzero 
anti-diagonal terms found in the experiment do not appear 
in the fit even though such terms may arise in our amplitude 
damping model. Due to the details of the fitting procedure, 
the optimal fit matches large elements closely but does not 
capture the detail of the small anti-diagonal. 
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FIG. 3. Determination of the strength of the scalar FIG. 5. Determination of the amplitude damping rates Tj . 

J-coupling between spins. The circles represent j = 1 (proton), and the squares repre- 

sent j = 2 (carbon). 




time (s) 



FIG. 4. Determination of the phase damping rates 7j . The 
circles represent j = 1 (proton), and the squares represent 
j — 2 (carbon). 73 is fitted to zero to within experimental 
uncertainty. 



J 220 ± 3 Hz 

7" 1 0.73 ± 0.03 s 

■y^ 1 0.21 ± 0.01s 

r^ 1 12.7±0.3s 

r^ 1 22.4 ± o.7s 



TABLE I. Dynamical time scales for CHCI3, as measured 
by QPT. Recall that the proton is spin 1 and the 13 C nucleus 
is spin 2. 73 is fitted to zero to within experimental uncer- 
tainty. 
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Our results are summarized in Table |. In addition 
to these data, we found the correlated phase damping 
rate 73 to be zero to within statistical precision. Unfor- 
tunately, insufficient precision was available to determine 
the amplitude damping temperature parameters fij . This 
is not a significant difficulty, as they are certain to be near 
i for high temperature systems. The data are consistent 
with this value. 

For comparison, the relevant decoherence parameters 
were measured by standard techniques. By the inversion- 
recovery technique, we found T± = 18.5 s for the proton 
and Ti = 21.1 s for carbon. These values agree roughly 
with the time scales T" 1 for amplitude damping. Us- 
ing the Carr-Purcell-Meiboom-Gill sequence, we found 
T2 = 4.7 s for the proton and T2 — 0.26 s for carbon. 
However, note that the relevant phase damping time scale 
is actually T£, the time scale for relaxation due to mag- 
netic field inhomogenity, which can be measured by fit- 
ting the free induction decay. By this method, we find 
T 2 * = 0.86 s for the proton and T 2 * = 0.20 s for carbon, 
values which agree more closely with the time scales 7" 1 
measured by QPT. 

As previously mentioned, the QPT preparation step is 
nonnegligible for the controlled-NOT gate fidelity experi- 
ment. However, this is not the case for the decoherence 
measurements: the time scales are long compared to the 
time to perform the QPT preparation, so small errors in 
the preparation are unimportant. For the experiments 
at short times to measure the parameter describing uni- 
tary evolution (J), the preparation is nonnegligible, but 
the operation is simple enough that we are still able to 
determine J to within two standard deviations. 



IV. CONCLUSIONS 

We have demonstrated the implementation of quan- 
tum process tomography to characterize the dynamics of 
a two-qubit NMR quantum computer. Such techniques 
should prove useful in the future for diagnosing quan- 
tum information processing devices. However, we should 
stress that due to the exponential size of %, QPT can 
only be used to characterize the dynamics of sufficiently 
small systems. For large systems, one may only be able 
to perform QPT on a small part of the total system, as- 
suming independence from the rest of the system — an 
assumption which can, in turn, be checked using QPT. 

Our measurements of gate fidelities underscore the dif- 
ficulty of implementing highly accurate logic gates in 
NMR. These fidelity measurements suggest a per-gate 
error rate of the order of 10 _1 , in agreement with previ- 
ous implementations of quantum algorithms |Hj|| . It has 
been suggested that error rates of 10 -2 can be achieved 
in solution NMR [p2[ . Indeed, if the previously dicussed 
cancellation of errors in a long sequence of gates leads to 
a T2-limited computation (as observed in |18|| ), then an 
optimistic estimate of T2 = 1 s and J = 100 Hz leads 



to an approximate error rate of 10 -2 . Nevertheless, even 
this falls far short of the requirements for fault-tolerant 
quantum computing Q : the least stringent estimates in- 
dicate a threshold for fault-tolerance of about 10~ 4 [M. 
Clearly, the development of fault-tolerant NMR quantum 
computers will require substantially modified techniques. 
One might construct an alternative model for fault toler- 
ance, perhaps based on topological properties [p5| . How- 
ever, with a conventional approach to fault tolerance, we 
will cither need significant improvement of gate fidelity 
over what is currently achievable, specialization of fault 
tolerant protocols to the specific errors which occur in 
NMR, or both. 

In turn, analyses of the threshold for fault-tolerant 
quantum computation typically involve an assumption 
of independent errors. Thus, for the experimental im- 
plementation of fault-tolerant computers, it is important 
to understand the error model in detail, and in partic- 
ular, the extent to which the independence assumption 
holds. The results presented in Section [II B show that 
this model is well- understood for the 13 CHCl3 NMR sys- 
tem, and furthermore, that the errors are uncorrelated. 
The data are fit well by a combination of phase damping 
and independent generalized amplitude damping, with 
the phase damping essentially uncorrelated. Thus there 
exist simple computational systems in which an uncorre- 
lated error model is reasonable. 

However, note that there are systems which do exhibit 
correlated decoherence, namely those that exhibit "cross- 
relaxation" (in the terminology of NMR). For example, 
the simple two-spin proton-carbon system provide d by 
sodium formate is known to possess cross-relaxation [ |26| . 
Investigation of such systems by QPT might be interest- 
ing not only for the purpose of quantum computing, but 
also for studying cross-relaxation in the context of con- 
ventional NMR. 
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